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REFLECTION  BY  DEFECTIVE  DIFFUSION  BONDS 


J.  D.  Achenbach  and  Y.  Xu 


Center  for  Quality  Engineering  and  Failure  Prevention 
Northwestern  University 
Evanston,  IL  60208-3020 


The  reflection  of  ultrasonic  wave  notion  by  planar 
defect  distributions  is  investigated  in  this 
paper,  with  a  view  towards  a  aethod  to  detect  and 
characterize  bond-plane  defects  and  with  an 
ultinate  objective  of  providing  information  on  the 
bond  strength  of  metal -to -metal  bonds.  The  paper 
is  primarily  concerned  with  an  approximate 
approach  to  analyze  reflection  by  cavities  and 
cracks  which  are  distributed  in  a  bond  plane.  The 
results  show  some  of  the  characteristic  features 
of  reflection  with  variation  of  such  parameters  as 
angle  of  Incidence,  frequency,  defect  size,  defect 
shape  and  defect  spacing. / 


I.  Introduction 

VMetal  to  metal  bonds  by  the  diffusion  welding 
process  are  becoming  more  frequently  used  for  the 
fabrication  of  complicated  parts,  such  as 
Integrally  bladed  compressor  and  turbine  rotors  in 
gas  turbine  engines.  Diffusion  bonds  may, 
however,  have  three  kinds  of  defects:  volumetric 
defects  (voids,  inclusions),  crack-like  defects, 
and  crystalline  grains.  ^ 

Defects,  when  they  do  occur,  are  generally 
distributed  in  the  plane  of  the  bond,  and  they 
tend  to  provide  a  plane  of  partial  reflection  for 
ultrasonic  wave  motion.  Reflected  and  transmitted 
waves  carry  Information  on  the  nature  of  the 
defects,  and  hence  the  ultrasonic  method  can,  in 
principle,  be  used  to  detect  and  characterize 
bond-plane  defects.  In  this  paper  reflection  and 
transmission  by  a  planar  distribution  of  one  kind 
of  defect,  namely  voids,  is  studied  in  some 
detail. 


Reflection  and  transmission  of  time-harmonic  waves 
by  planar  arrays  of  cracks,  cavities  and 
inclusions  have  been  studied  by  Achenbach  et  al, 
[l]-(3],  [7]  and  Thompson  et  al,  [5] -[6].  Angel 
and  Achenbach,  [1],  and  Hikata  and  Achenbach,  [2], 
have  obtained  numerically  rigorous  reflection  and 
transmission  coefficients  for  the  two-dimensional 
problem  of  a  periodic  distribution  of  cracks. 
Nonperiodic  crack  distributions  have  been 
considered  by  Sotlropoulos  and  Achenbach  [3] -[A], 
Thompson  and  co-authors,  ( 5 J - (6 J ,  have  studied 
reflection  from  imperfect  bonds,  both  for  one- 
dimensional  and  tvo-dimenslonal  planar  crack 
distributions  based  on  the  use  of  a  distributed 
spring  model.  Reflection  and  transmission  by  a 
periodic  array  of  spherical  cavities  has  been 
Investigated  by  Achenbach  and  Kitahara  [7], 
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The  relevant  parameters  for  the  present  problem 
are  the  angle  of  incidence,  the  frequency,  length 
measures  of  the  spacing  of  the  cavities  and  their 
sizes,  and  the  properties  of  the  material. 

Specific  results  are  presented  for  reflection  by  a 
distribution  of  spherical  cavities.  For  the 
special  case  of  a  periodic  distribution  the 
reflection  and  transmission  coefficients  computed 
as  functions  of  the  frequency  agree  very  well  vith 
the  exact  numerical  results  of  Ref.  [7]. 


Fig.  1  Distribution  of  cavities  in  plane  "  0 

ii-  Breblsa  ForaulaUan 


The  plane  defined  by  x.  -  0  in  an  unbounded, 
homogeneous,  isotropic;  linearly  elastic  solid 
Intersects  a  distribution  of  cavities  of  arbitrary 
shape.  There  are  no  cavities  that  are  not 
Intersected  by  the  plane  x,  -  0.  For  the  present 
discussion  we  consider  an  area  of  the  plane  x,  -  0 
defined  by  |xj|  £  1,  |x-|  v.  This  area 
intersects  M  cavities,  which  are  numbered  m  -  1, 

2  •••  M.  The  geometry  is  shown  in  Fig.  1. 


A  plane  longitudinal  time-harmonic  wave  is 
incident  on  the  distribution  of  cavities.  The 
time-harmonic  term  exp(-lut),  where  u  is  the 
angular  frequency,  will  however  be  omitted.  The 
Incident  longitudinal  wave  is  of  the  general  form 
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Ad^exp(ikjjjL*x) 


Here  A  is  the  amplitude,  and 


.L  L 
d  -  g 


(0,sin#L,cos*L) 


(1) 


(2) 


□ 

i) 


where  I.  is  the  angle  of  incidence.  The  surfaces 
of  the  cavities  are  free  of  tractions.  Thus  for 
the  m-th  cavitv  we  have 
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x  <  S^:  t^  -  “  0  .  (3) 

where  r. ,  denotes  Che  stress  censor,  end  n  Is  the 
unit  vector  normal  to  Che  cavity  surface  3^. 

It  Is  assumed  that  far  from  the  plane  Intersecting 
the  cavities  (Xj  -  0) ,  the  dominant  vave  fields 
may  be  written  as  reflected  and  transmitted  plane 
longitudinal  and  transverse  waves  of  the  forms : 

reflected:  u®0  -  RaA  d80  exp( (4) 

transmitted:  u7<*  -  T0 A  d7°  exp(ikog7a*x)  (5) 

In  R0,  we  have  a  -  L,  TV  and  TH  for  reflected 
longitudinal,  reflected  vertically  polarized 
transverse  waves,  and  reflected  horizontally 
polarized  transverse  waves,  respectively. 

Analogous  definitions  hold  for  T0,  vhere 
a  -  L,  TV  and  TH.  The  unit  vectors  appearing  in 
Equation  (4)  are 

g^WO, sin*L,  -cos*L) ,  E«V“-(0,si»#T.  'co* V 
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where  a  -  L,TV,TH,  i  and  w  are  the  length  and  the 
width  of  a  rectangular  area  In  the  plane  x.  -  0, 
and  M  is  the  number  of  cavity  intersections  In 
that  area.  Also 
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for 

the  m-th 

the 

cavity  volume 

(6a, b) 

dRL  RL^ dRTH_(1  > 0  0)  _  ,  Cos<_ . sln«_) 

1  1(7a,b,c) 

The  angles  #L  and  I T  are  related  by  Snell's  law: 

^•In^L  -  Vin*T  <8) 

The  unit  vectors  in  Equation  (5)  are  defined  as 

£*VV,  (0.sin«T.cos#T)  (9a, b) 

d^-d.O.O),  d777  -  (0.  -cos»T,sln»T)  (10a, b) 

Far  from  x-  -  0  the  total  displacement  fields  may 
now  be  written  as: 

x3  >  0:  u(x)  -  J  T®  A  d70exp(lkag7“«x) , 

a  -  L.TV.TH  (11) 

x3  <  0:  u(x)  -  u1  +  E  R°  A  d*0exp(lkag8a«x) , 

a  -  L.TV.TH  (12) 

Equations  (11)  and  (12)  Imply  that  the  plane 
containing  the  cavities  is  taken  as  a  plane  of 
homogeneous  reflection  and  transmission.  It  has 
been  shown  by  Achenbach  and  Kitahara  [7]  that  this 
is  rigorously  correct  for  the  lowest  modes  when 
the  cavities  are  spheres  and  are  periodically 
distributed  in  the  plane  x-  -  0.  For  an  arbitrary 
distribution  of  cavities  Equations  (11)  and  (12) 
are  assumptions  which  should  be  good  at  low 
frequencies  and  at  some  distance  from  the  plane 


V^j*0  -  f  Ujttj  exp(-ikQgR0*x)dS 


(17) 


_  m 

with  an  analogous  definition  for  V^j  ,  while 
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hi.  Cavity  Interaction  Factor 


(18) 


*ij  ■  Pi“dja  +  Pj“di“  for  0  -  TV-™.  (19) 


Equations  (13) -(16)  show  that  the  reflection  and 
transmission  coefficients  for  a  planar 
distribution  of  cavities  can  be  written  in  terms 
of  the  cavity  volume  tensors  V^j70  vhere  7  -  R 

(reflection)  or  7  -  T  (transmission)  and 
a  -  L,TV  or  TH.  An  approximate  calculation  of  the 
reflection  and  transmission  coefficients  can  be 
based  on  the  use  of  a  cavity  interaction  factor  to 
calculate  the  cavity  volume  tensors.  Let  us  first 

consider  the  effect  on  reflection  and  transmission 
coefficients  of  the  interaction  between  two 
neighboring  cavities.  These  cavities  are  denoted 
by  1  and  2. 


For  cavity  1,  the  displacement  components  u!  can 
be  written  as  1 


u‘  -  u«°  +  u‘* 


(20) 


As  shown  by  Achenbach  and  Kitahara  [7]  and 
subsequently  in  more  general  form  by  Sotiropoulos 
and  Achenbach  (3],  the  reflection  coefficients  R 
and  the  transmission  coefficients  T°,  can  be 
obtained  by  application  of  the  reciprocal 

identity.  For  details  we  refer  to  [3]  and  [7j. 
The  results  are 


vhere  ui 0  are  the  displacement  components  induced 
by  the  Scattered  field  in  the  absence  of  the 
second  cavity,  and  u£J  are  the  displacement 
components  on  cavity1!  due  to  the  presence  of 
cavity  2.  The  analogous  components  of  the  cavity 
volume  tensor  can  be  written  as 

v10,7<>  v12,7<> 
ij  ij 


V1’70  - 
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(21) 
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defined  by  Equation  (17),  and  V*°'7®  and  V j j • 70 

12  TO 

are  defined  analogously.  To  determine  V  ,T  It 
Is  convenient  to  define  the  cavity  interaction 
factor 

^j’7®  by  the  relation 

<»> 

12  to 

Th«  cavity  Interaction  factor  gives  the 

ljotn 

components  of  the  cavity  volume  tensor  of  cavity  1 
due  to  the  presence  of  cavity  2,  in  terns  of  the 
components  of  the  total  cavity  volume  tensor  of 
cavity  2. 
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Details  of  the  calculations  of  ,7^  and 

flj?’7®  are  given  by  Xu  and  Achenbach  [8].  The 
lj  nn 

12  ufl 

Interaction  factor  will  depend  on  a  number 

lj  nzi 

of  parameters.  Prominent  among  these  are  of 
course  the  shape  of  cavity  1,  the  shape  of  cavity 
2,  the  distance  between  the  cavities  as  compared 
to  their  dimensions,  the  angle  of  Incidence,  and 
the  frequency  of  the  Incident  wave. 

The  concept  of  the  Interaction  factor  can  now 
easily  be  extended  to  a  distribution  of  cavities, 

l.e. ,  for  cavities  k  and  s  we  define  by  the 

relation 

'C.1"  <«> 

where  the  summation  Is  over  m  and  n  only.  Then, 

with  the  additional  assumption  that  -Sr  Is  not 

affected  by  the  presence  of  cavities  other  than  k 
and  s,  we  can  write 

.  I  <“> 

lJ  s-l(s^k) 

Once  the  cavity  Interaction  factors  are  lcn  a 
relatively  simple  way  of  calculating  the 
components  of  the  cavity  volume  tensor  is  by 
use  of  an  Iteration  procedure. 

For  cavity  k  the  components  of  V^j7®  are  given  by 

Equation  (24).  In  first  approximation  V^7®  Is 

taken  as  V^'7®,  l.e.,  the  value  of  V^7®  in  the 

absence  of  all  the  other  cavities.  The  value  of 
V^j7®  so  obtained  Is  termed  V^j’7  ,  l.e., 

vkl.7«*  _  yko.T*  +  j  d^’7®  V*0’7®  (25) 

U  s-l(s.k)  « 
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ySL,70 


s-l(sHk) 


In  the  next  step  V*j’7®  1*  calculated  by  Equation 

(25),  and  the  result  is  substituted  In  the  right- 
hand  side  to  vield 


This  process  can  be  repeated  as  many  times  as 
necessary.  The  procedure  Is  particularly  useful 
when  the  array  of  cavities  la  a  periodic  one, 
because  then 

V^7®  -  vj}*7®  (27) 

Is  rigorously  correct. 


The  simplest  result  is  achieved  when  the 
Interaction  between  neighboring  cavities  Is 

neglected.  For  that  case  the  summation  of  terms 
on  the  right-hand  side  of  Eq.(24)  disappears.  The 
expression  for  the  coefficient  of  reflection, 
given  by  Eq.(13)  reduces  to 


„o  _  _ M _  yO.fc* 

R  2coafa  *1J  Vij 


Here  M  -  H/tv  denotes  the  number  of  cavity 
Intersections  per  unit  area  In  the  plane  x^  -  0, 
and 

g 

v°,Ra  _  1  r  v°,R®  (29) 

Vij  M^iVlj 

are  the  components  of  the  averaged  cavity  colume 
tensor. 

For  a  distribution  of  spherical-cavities  of  radius 
d,  and  for  normal  Incidence,  |R  |  according  to 
Eq.  (13)  has  been  plotted  versus  lLd  In  Fig.  2. 

Here  the  Poisson's  ratio  of  the  material  was 


Here  the  Poisson’s  ra 
selected  as  u  -  0.25. 
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Fig.  2  Reflection  coefficients  for  spherical 
cavities 

It  is  of  interest  to  compare  |RL|  for  the  case  of 
no  cavity  Interaction  with  corresponding  results 
for  the  case  where  Interaction  has  been  taken  Into 
account.  Achenbach  and  Kitahara  [7]  have 
presented  rigorous  results  for  reflection  by  a 
periodic  array  of  spherical  cavities  whose  centers 
are  spaced  distances  a  and  b  in  the  x.  and  x_ 
directions,  respectively.  For  a/d  -  5/d  -  3,  [R  | 
according  to  Ref. [7]  has  also  been  plotted  in  Fig. 
2.  It  is  noted  that  cavity  interaction  starts  t” 
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as  Che  frequency  becomes  larger.  Also  plotted  In 
Fig-  2  is  the  reflection  coefficient  for  cavity 
interaction,  but  where  the  components  of  Che 
cavity  volume  tensor  have  been  computed  according 
to  the  Iteration  procedure  of  Eqs. (25)-(27) . 


Fig.  3  Oblique  incidence  on  periodic  distribu¬ 
tion  of  spherical  cavities 

Figure  3  compares  the  reflection  and  transmission 
coefficients  for  an  angle  of  incidence  t.  -  30° 
computed  by  the  method  of  this  paper  with  the 
results  of  Ref. [7].  Good  agreement  is  again 
observed. 


V.  Statistical  Distribution  of  Cavities 


Next  we  consider  the  effect  on  a  central  cavity  of 
a  random  distribution  of  neighboring  cavities. 

Let  the  neighboring  cavities  that  have  an  effect 

on  the  central  cavity  be  distributed  in  an  annulus 
of  Area  A,  defined  by  <  r  <  a  around  the 

central  cavity,  see  Fig.  4.  The  probability 
density  function  of  the  neighboring  cavities  is 
defined  as 


P<r.*) 


(30) 


l.e.,  the  distribution  is  completely  random  in 
both  the  radial  and  circumferential  directions. 


Fig.  4  Configuration  of  central  cavity  and 
neighboring  cavity 


in  the  usual  manner.  For  example, 

.ks.ja  _  1  2*  *- 

^ijmn  2s  (a 


‘min 


)  i  r 

“ln  (31) 

The  cavity  interaction  factor  naturally 

depends  on  the  distance  between  the  centers  of  the 
cavities.  After  an  XjXj-aystem  of  axes  to  which 
the  Indices  l.j.rn  and^n  refer,  has  been  chosen, 

^ijmn*  wil1  M  *lso  'kP81*1  on  che  angle  +  between 

the  line  connecting  the  centers  of  two  cavities 
and  the  Xj-axis,  see  Fig.  4.  The  relation  between 

^ijmn*  (r>^)  *nd  (r,^)  is,  however,  of  the 

following  simple  form: 


^  <*•*>  “  QsiW/^V.WtjW 


where 


[QJ 
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c  osj 
sin^ 
0 


sin^ 

cos 

0 


(32) 


(33) 


Substitution  of  Eq.(32)  into  (31)  allows  a  simple 

evaluation  of  the  components  of  <  >‘ln 

ij  am 

terms  of  (r,0). 

In  the  next  step  of  the  argument  it  is  realized 
that  for  a  totally  random  distribution,  every 
cavity  can  be  considered  as  the  central  cavity, 
l.e.,  all  cavities  are  the  same.  An  application 
of  the  Integral  of  Eq.(31)  to  Eq.(24)  then  leads 
to  the  relation 

Vjj 70  -  ^j0’70  +  <«-!>  <  (34) 

Here  it  has  been  taken  into  account  that 

an 

*nd  Vij  *re  (■ndeP«"dent  of  r  and  Also,  H  is 

the  total  number  of  cavities  in  the  area  wa2 

max' 

Substitution  of  Eq. (34)  into  Eqs. (13)  -  (16) 
yields  the  corresponding  expressions  for  the 
reflection  and  transmission  coefficients. 

Numerical  calculations  have  been  carried  out  for 
specific  choices  of  a 


min’ 


define  the  quantity  e  as 

«  -  -i&i 

*  a2 


max 


and  M.  Ve  also 


(35) 


i.e.,  c  defines  the  ratio  of  the  area  of  the 
cavity  cross-sections  to  the  total  area  that 
Influences  the  central  cavity. 


For  Poisson's  ratio  v  -  0.25,  a  ,  /d  -  2.8  and 

Bln 

c  -  0.126,  Fig.  5  shows  the  reflection  and 
transmission  coefficients,  for  the  case  of  normal 
incidence. 
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Reflection  and  transmission  coefficients 
for  random  distribution  of  spherical 
cavities 


The  analysis  of  this  paper  Is  also  valid  for 
crack- like  flavs  in  the  plane  x,  -  0.  For  crack- 

«  ••  *—%  _ «  _  _  _  a  r _ i  _  _  i  *nv  — 


like  flaws  we  have 


"l  “  n2  = 


—  0  and 


expressions  for  the  reflection  and  transmission 
coefficients  become  simpler  In  form,  since 


V11  "  V22 


-  V31  -  V32  =  0 


The  remaining  components  of  the  cavity  volume 
tensor  now  become  components  of  the  crack- opening 
tensor.  Further  details  have  been  worked  out  In 
Ref. [4] ,  where  It  has  also  been  shown  that  the 
Mode-I  stress -intensity  factor  for  a  distribution 
of  equal -sized  cracks  can  be  directly  related  to 
reflection  data. 


This  work  was  carried  out  in  the  course  of 
research  sponsored  by  the  Office  of  Naval  Research 
under  Grant  No.  N00014-89-J-1362  to  Northwestern 
University. 
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